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1 Introduction 

We consider a heterogeneous porous medium supplied by several reservoirs of a fluid, 
represented by a bounded domain R of M” with locally Lipschitz boundary dR = 
Fi U r 2 , where Fi is the impervious part of the boundary, F 2 is the part in contact 
with either air or the fluid reservoirs. 
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The fluid infiltrates through hi obeying to Darcy’s law 


V = -a(x).V(u + Xn), 

where a{x) = {aij{x))ij is the n X n permeability matrix of the medium, x = 
(xi, ...,Xn), V is the fluid velocity and u its pressure. 

We are concerned with the problem of finding the pressure u and the saturation x 
of the fluid inside D. Using the mass conservation law, Darcy’s law, the boundary 
conditions and the initial data, we obtain the following strong formulation for our 
problem (see 131): 


^r>0, 0<x< 1, u(l-x) = 0 

in Q 

div{a{x){'Vu + xe)) — {au + x)t = 0 

in Q 

u = (p 

on S2 

{au+ x)i-,0) =auo + Xo 

in D 

a{x){Vu + xe) • = 0 

on Si 

a(x)(Vn + xe) ■ v < 0 

on S4 


where a, T are positive numbers, Q = Q x (0, T), Si = Ti x (0, T) is the impervious 
part of 5 D, S2 = r2 X (0,T) is the pervious part, S3 = S2 H {()) > 0} is the part 
covered by fluid, and S4 = S4 n {(/> = 0} is the part where the fluid flows outside D. 
0 is a nonnegative Lipschitz continuous function defined in Q, is the outward unit 
normal vector to 9D, e= (0,...,0,1) € M”, uq,xo ■ ^ are functions satisfying 

for a positive constant M 

0 < uq{x) < M, 0 < Xo(2:) < 1 for a.e. x € D. (1.1) 

For a(x), we assume that we have for two positive constants A and A 

€ M"", for a.e. x G Q < o(x).^.^. (1.2) 

G M”, for a.e. x G fl |a(x).^| < A|^|. (1.3) 
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Moreover, we assume that 


div{a{x)e) G L^(Q). 


(1.4) 


Using the strong formulation, we are led to the following weak formulation 


Find (u,x) € L^(0,T; x L°°{Q) such that : 

(i) u > 0, 0 < X < 1; ^^-(1 — x) = 0 a.e. in Q 

(ii) u = (j) on S 2 

(in) / [a(x) (Vtt + xe) • — (au + x)^i] dx (it 

Jq 

< / (Xo+ a^io)?(3:,0)(ix 

Jn 

G H^iQ), C = 0 on S 3 , ^ > 0 on S 2 , 

^{x,T) = 0 for a.e. x G 0. 


For the existence of a solution of the problem (P) in the homogeneous case (a(x) = 
In), we refer to [2] and [T] respectively in the incompressible (a = 0 ) and compress¬ 
ible (a > 0) cases. For the heterogeneous case, we refer to |16) in a more general 
framework under assumptions (1.1)-(1.4) for both incompressible and compressible 
cases. For the incompressible case with nonlinear Darcy’s law, we refer to m, m 
and [3] respectively for Dirichlet, Neuman and generalized boundary conditions. Re¬ 
garding regularity of the solution, we refer to [ 2 ] and [T], where it has been proved 
when a(x) = In that y G C^([0,T]] LP(Il)) for all p > 1 in both incompressible and 
compressible cases, and that u G (^^([O,T];LP(D)) for all 1 < p < 2 in the com¬ 
pressible case. Extensions to the quasilinear and incompressible case were obtained 
in HU, m and m in both homogeneous and nonhomogeneous frameworks. The 
authors of this paper were recently able to extend the above regularity result in |14| 
to a more general framework under weaker assumptions on the data. 
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In this paper, we are mainly concerned with the uniqueness of the solution of the 
problem (P). This question was first addressed for a rectangular homogeneous dam 
in m and [7] respectively for a formulation based on quasi-variational inequalities 
and for the formulation (P) with a dam wet at the bottom and dry near to the top 
in the second case. Uniqueness of the solution for a homogeneous dam with general 
geometry was established by the method of doubling variables in [T], but it is not 
obvious wether it works in the heterogeneous situation. Extensions to a quasilinear 
operator modeling incompressible fluid flow governed by a nonlinear Darcy’s law with 
Dirichlet, or Neuman boundary conditions were obtained in im and m respectively. 

Our main result in this work is the uniqueness of problem (P) solution for a heteroge¬ 
neous and rectangular porous medium assuming it remains wet near the bottom and 
dry close to the top. Our method is inspired by an idea from [7] in the homogeneous 
case and relies on solution regularity that has been recently obtained in M- Our 
uniqueness result is new in the heterogeneous and rectangular framework, but most 
likely the technic is limited to that particular geometry like shape. 

2 Preliminary Results 

In this work we shall be interested in the following situation of a two dimensional 
rectangular dam D = (0, L) x {0,K), with L,K >0, and Ti = [0, P] x {0}, r 2 = 
({0} X [0,K]) U ([0,P] X {K}) U ({P} X [0,P:]) (see Figure 1). 

We also assume that 


a{x)e G 

(2.1) 

div{a{x)e) > 0 a.e. in D, 

(2.2) 

00 < 0 < 01 on S 2 , 

(2.3) 


where <J)q and cpi are two nonnegative Lipschitz continuous functions dehned on D 
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Figure 1 


and satisfying for some eo > 0 small enough 

' 

'/>o( 0 , X2) = (po{L, X2) = (eo - X2)'^ 

' 01 (0, X 2 ) = (j)i{L, X 2 ) = {K -eo- X2)+ (2-4) 

(l)o{xi,K) = (l)i{xi,K) = 0 . 

Let us now denote by [vi,'yi) the solution of the stationary problem corresponding 
to 0 j, i = 0,1 (see |5]) 


(Pf) 


Find (ui, 7 i) € x L°°{Q) such that : 

(z) Uj > 0 , 0 < 7 j < 1 , Uj.(l — 7 j) = 0 a.e. in n 

(m) Vi = (t)i on r 2 
(in) / a(x)(Vuj + 7 je).V.^dx < 0 

Jn 

^ ^ = 0 on r 2 n { 0 i > 0 }, ^>o on r 2 n{ 0 j = o}. 


We recall the following uniqueness result. 

Theorem 2.1. Assume that (2.1)-(2.2) hold. Then the solution {vi,'ji) of (Pf) is 
unique and satisfies 

7i = X{vi>0}- (2.5) 
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Proof. When a{x)e € (7^(0), we refer to [1], since we cannot have pools in a rectan¬ 
gular dam. One can also argue as in m- When a(x)e G one may combine 

Theorem 5.1 of [5] and the proof of Theorem 6.3 in |4] to establish the result. 

□ 


The following properties for the solution of {Pf) (i = 0,1) hold. 
Theorem 2.2. Assume that (2.1)-(2.4) hold. Then we have 


(eo - X2] 

)^ < 'Oo < < {K — eo — X2)^ a.e. in 

(2.6) 

70 = 1 

o.e. in fl n {0 < 3:2 < eo} 

(2.7) 

7i = 0 

a.e. in Q Cl {K — eo < X2 < K}. 

(2.8) 


Proof. First, we remark that for each k G (0,K), (fe — X 2 )"*~ satisfies the equation 

/ a{x){V{k - X 2 )^+ X{{k-x 2 )+>o}e)-'^Cdx = a{x){-e + e).V(,dx = 0. 

Jo. JQ.r\{x2<k} 

It follows that {{k — * 2 )"*", X{(A:-a; 2 )+> 0 }) is a solution of the stationary dam problem 

for the boundary Dirichlet data {k — X 2 )'*' on r 2 . Adapting the proof of uniqueness 

in [1] and arguing as in |10) . and using (2.1)-(2.4), we obtain 


(eo - X2)'^ < vq < vi < {K - eQ - X2)'^ in 
X{{eo-x2)+>o} < 70 < 7i < X{{K-eo-x2)+>o} a.e. in Q. 


Hence (2.6)-(2.8) follow. 


□ 


Remark 2.1. Theorems 2.1 and 2.2 remain true without the regularity assumption 
(2.1) provided the following assumptions on the permeability matrix hold (see fWf) 


012 = 0 
da22 


dxo 


> 0 


a.e. in Q, 
in T>'{Q,). 
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Next we will construct a solution corresponding to a dam that is wet up to X2 = cq 
and dry above X2 = K — eq over the whole interval [0, T]. 

Lemma 2.1. Assume that ( 2 . 1 )-( 2 . 4 ) hold and the initial data satisfies 

vq < uq < vi a.e. in (2.9) 

7o < Xo < 7i tJ-e. in (2.10) 

Then there exists a solution {u,x) of problem (P) such that 

vq <u <vi a.e. in Q (2-11) 

7o < X < 7i inQ. (2.12) 


Proof. Let Vi^ be the solution of the approximating problem of the stationary problem 


(Pi), * = 0,1 


(PL 


Find Vif_ G H^{Q) such that : 

(i) Vie = (t^i on r2 

{ii) / a{x){Vvie + He{vie)e).'Vfdx = 0 

Jn 


e = 0 onFs, 

\ 


where Hfis) = min(l, s'’'/e) is an approximation of the Heaviside graph H{s) = 
[ 0 , l]X{o} + X(o ,oo) * 

Let Ue be the solution of the following approximating problem of the problem (P) 


(P.) 




Find Ue G H^{Q) such that : 

(i) Ue = 4> on T ,2 

{a) / [a{x) (Vue + iLe('Ue)e) • + eUetit - Ge{ue)ii\ dxdt 

JQ 


+ / Ge{ue{x,T))^{x,T)dx = / {auQe{x)+ Xo^{x))i{x,{))dxdy 
J Q, J Q 

y^eH\Q), e = 0on S2, 
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where = min(uo,fie) and xoe = ram{xo, H^{vu))- 
If ^ G H^{Q), ^ = 0 on S 2 , we have from {PfJ{ii) : 


[ [a{x) {Vvie + H^{vie)e) ■ - (avie + H^{vie))^t] dxdt 

JQ 

+ / {avu +H^{vie))i{x,T)dx = / {avu + H^{vu))i{x,Q)dx. 

Jn Jn 


(2.13) 


(^U V 

For h > 0, the function ^ ^- belongs to H^(Q) and satishes = 0 

Ue - Vie 

on S 2 since (j) < 4>i on S 2 . Writing (2.13) and {Pe){ii) for ^ and subtracting the 
two identities from each other, we get by taking into account (2.9)-(2.10) 


IQ 


l(x)(V(Ue - Vie) + {HeiUe) - He{vie))e) .V+ ^iPe - Vle)t^St 
-{a{ue - Vie) + He{ue) “ Htivie))^st dxdydt 
+ [ {a{Ueix,T) - Vie) + He{Ue{x,T)) - He{vie))^six,T)dx 

Jn 

= [ {a{u 0 e - Vie) + XOe - He{vie))^5{x,0)dx <0. (2.14) 

Jn 

By Lemma 2.1 of US], we obtain from (2.14) 


Ue < Vie a.e. in Q 


(2.15) 


and by the monotonicity of He, we get 


He{ue) < Heivie) a.e. in Q. (2.16) 

We recall that from the proof of existence (see [8] or |16| for example), we know that 
we have up to a subsequence 

Ue ^ u weakly in L^(0, T; Lf^(n)). (2-17) 

He{ue) X weakly in L?{Q) (2.18) 


where {u,x) is a solution of problem (P). 





Similarly, we have since the solution of problem (Pf) is unique 


Vie vi 
He{vie) 7l 


weakly in 
weakly in 


(2.19) 

( 2 . 20 ) 


Now, let ^ € P(Q) with ^ > 0. Passing to the limit, we obtain by using (2.15)-(2.20) 
/ (ui — u)^dxdt = lim / (vie — Ue)^dxdt > 0, 

Jq 

[ (71 - x)idxdt = lini [ - He{ue))^dxdt > 0, 

Jq 

which leads to 


u < vi a.e. in Q 
X < 7 i a.e. in Q. 
(uoe -u^- 5)+ 


( 2 . 21 ) 

( 2 . 22 ) 

belongs to (Q) and satisfies 


Similarly, for <5 > 0 the function ^,5 = 

t>0e - Ue 

= 0 on S 2 since 4>q < (p on S 2 . Then by taking into account (2.9)-(2.10), we get 


'Q 


l{x) (V(uoe - Ue) + {HeivQe) “ He{Ue))e) .V^S + e{voe “ Ue)tCst 
-{a{v 0 e - Ue) + Heivoe) “ He{Ue))Cst dxdt 
+ [ {a{v 0 e - Ue{x,T)) + H^{v 0 e) - He{Ue{x,T)))^six,T)dx 

Jn 

= f {a{voe - Uoe) + H^{voe) “ X 0 e)^s{x, 0)dx < 0. (2.23) 

Jo. 

By Lemma 2.1 of [16], we obtain from (2.23) 


< Ve a.e. in Q 


(2.24) 


and by the monotonicity of He, we get 


He{voe) < H^{ue) a.e. in Q. 


(2.25) 
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Arguing as above and using (2.24)-(2.25), we obtain by passing to the limit up to a 
subsequence, that we have for any ^ G 'D{Q) with ^ > 0 

(u — vo)^dxdt = lim / (u^ — VQf)^dxdt > 0, 

JQ 

[ (x -'yo)^dxdt = lim f - H^{voe))Cdxdt > 0, 

Jq 

which leads to 


vq < u a.e. in Q (2.26) 

7 o < X a.e. in Q. (2.27) 

Combining (2.21)-(2.22) and (2.26)-(2.27), we obtain (2.11)-(2.12). □ 

Remark 2.2. Assume that a{x)e € C''^’^(n). Then we get from (2.11)-(2.12) taking 
into aceount (2.5) 

u{x,t) >0 if 0 < X 2 < eo (2.28) 

u{x,t) = x{x,t) = 0 if K — eQ < X2 < K. (2.29) 


3 Uniqueness of the solution in rectangular dams 

In this section we assume that 

a(x) G C°’^(0), with iV = sup loo- (3.1) 

ij,k 

a{x) is a symmetric matrix. (3.2) 

Here is our main result. 

Theorem 3.1. Assume that (2.2) and (3.1)-(3.2) hold. Then the solution of the 
problem (P) associated with the initial data {uq,xo) o,nd satisfying (2.28)-(2.29) is 
unique. 
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Let s-iid {u 2 ,X 2 ) be two solutions of the problem (P) satisfying (2.28)-(2.29). 

Set 

w = ui — U 2 and r] = aw + Xi ~ X 2 - 


We consider the following problem 


Find n G -L^(0, T; Lr^(n)) such that : 
div{a[x)Vv) = —r] in Q. for each t G [0,T] 
V = Q on r2 
a{x)'Vv.i2 = 0 on Pi. 

The we have 


(3.3) 

(3.4) 

(3.5) 


Lemma 3.1. There exists a unique weak solution of the problem (3.3)-(3.5). 


Proof. First, we observe (see [H]) that aui + x* G C^{[0,T]; i = 0,1. As a 

consequence, we have r] G C^{[0,T]-, L'^{Q)). Let V = {v ^ H^{VL) / u = 0 on r 2 }. 
Then IL is a Hilbert space under the norm, and by applying Lax-Milgram’s 

Theorem and taking into account (1.2)-(1.3), there exists for each t G [0,T] a unique 
solution v{x,t) of the following problem 

' vi.,t)eV 

< ^ ^ (3.6) 

/ a{x)\/v{x,t).Vf,dx = / r]{x,t)f,dx G H. 

Jn Jn 

Choosing G P(0) in (3.6), we obtain (3.3) in 'D'{Q.) and therefore in C''^([0, T]; iF'iPt)) 
(3.4) is satisfied in the trace sense and (3.5). Writing (3.6) for ^ G C'°°(n) with = 0 
on r 2 , and taking into account (3.3), we obtain (3.4) in 

Choosing u as a test function in (3.6) and using (1.2), Holder and Poincare’s inequal¬ 
ities, we obtain 


|V?;(x,t)P(ix <[ \r](x,t)\‘^dx. 

Jn 


(3.7) 
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Integrating (3.7) over the interval [0,T] and using the fact that Xi ^ > Ui G 

L°°{0,T] L°°{Q)) (see [12)) we obtain 

[ \Vv{x,t)\‘^dxdt < ^ [ [ \r]{x,t)\‘^dxdt — (o,r,L (n)) ^ 

Jq Jo Jn ^ 

Using Poincare’s inequality, we obtain v € L^(0,T; Hence v is the unique 

solution of (3.3)-(3.5). 

□ 

Remark 3.1. By the regularity theory (see m for example), the solution v of the 
problem (3.3)-(3.5) satisfies v G L^(0, T; (^^(Huf lUf 2 )) nL^(0, T; R^(Huf 1 uf 2 )). 


Now, let us denote by g the mean with respect to t of a function g(x, t) defined by 

rt+h 

9 {x,t) = - J g{x,s)ds. 


Then we have 


g ^ g as h —>■ 0 

^ = ^ig{x,t + h) - g{x,t)). (3.8) 

Moreover, it is easy to check that 

div{a{x)Vv) = —rj in H for all t G [0, T] (3.9) 

V = 0 on r 2 (3.10) 

a{x)Vv.u = 0 on Pi. (3.11) 


Since y* = 1 (i = 1, 2) in a neighborhood of Pi, we obtain from {P)iii) 

{aui + Xi)t = div{a{x){Vui + y^e)) in V'{Q), 
a{x) {y Ui + e) m = 0 on Pi, z = l,2. 
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Using the fact that ui = U 2 on S 2 , and writing the previous two equations for (ui, xi) 
and {u 2 ,X 2 ) and subtracting them from each other, we get 


{aw+ xi - X2)t = div{a{x){Vw + {xi - X2)e)) mV'{Q). (3.12) 


w = 0 on S 2 
a{x){'\/w).v = 0 on Si. 


(3.13) 

(3.14) 


Then we have 


Lemma 3.2. For h small enough we have 


rjt = div{a{x){Vw + {xi - X 2 )e)) in V'{Q) 
w = Q on Ti 2 
a{x)Vw.v = 0 on Si. 


(3.15) 

(3.16) 

(3.17) 


Proof. (3.16) and (3.17) are a direct consequence of (3.13) and (3.14). To establish 

(3.15), let C G PiQ) such that for some tq > 0, supp{C) C H x (to,T — tq). We 

- 1 /■* 

denote by f the function defined by C(®)0 ~T C{x,s)ds. Since for \h\ < ro/2, 
the functions ±(^ are test functions for problem (P), we obtain for i = 1, 2 

/ a{x){Vui + Xi^) ■'^Cdxdt = / {aui + Xi)Ctdxdt. (3.18) 

JQ JQ 

For the right hand side of (3.18), we have by using change of variables 


'Q 


{aui + Xi)Ctdxdt = f {aui + Xi)TC{x^t)dxdt — [ {aui + Xi)-rC{x,t — h)dxdt 


'Q 


h 


'Q 


h 


= / {aui + Xi){x,t)\c,{x,t)dxdt - f [ ^{aui{x,t + h) + Xi{x,t + h))C{x,t)dxdt 

JQ h Jq h 

rT 


= {aui + Xi){x,t)YC{x,t)dxdt - [ [ + h) + Xi{x,t + h))C{x,t)dxdt 

JQ h Jn Jo h 


'Q 


[{aui{x,t + h) + Xi{x,t + h)) - {aui + Xi){x-,t)]C,{x-,t)dxdt 


= - {aUi + Xi)tCix,t)dxdt. 
JQ 


(3.19) 


13 


For the left hand side of (3.18), we have by integrating by parts 

/ a{x)(yui +Xie) ■'^(dxdt = - / a{x) [VUi + Xie) ds 

JQ Jq ^ Jo 


SI Qtdxdt 


Q '-Jo 


a(x) (Vui + XiG)ds ■ V[ — {C{x, t) — ({x,t — h))'^dxdt 


/ a{x) (siU i + XiG-)ds ■ SIC{x,t)dxdt 

Jo 

+ / — f a(x) (siU i + Xis) ds ■'VC(x,t — h)dxdt 

Jq d I Jo 

— I T I a(x) (Slui + Xi^)ds ■ SIC(x,t)dxdt 

Jq d I Jo 

f 1 r 1 

+ -r a(x)(SIui + Xie)ds ■ SIC(x,t)dxdt 

JnJo dlJo 

f 1 r 1 

/ — / a(x){yui-\-Xie)ds ■V((x,t)dxdt 

Jq d\-Jt -I 

/ a(x){yui-\-Xie) ■SJ(l(x,t)dxdt. 

Jq 


(3.20) 


Then we deduce from (3.18)-(3.20) that 

/ a(x)(SIui-\-Xic) ■ SIC(x,t)dxdt = — / (aui + Xi)tC(x,t)dxdt. 

Jq Jq 

Writing the last equation for i = 1, 2 and subtracting the two equations, we get 
/ a(x)(SIw + (Ja■ SIC(x,t)dxdt = — rjtC(x,t)dxdt 

Jq Jq 

which is (3.15). □ 

To prove Theorem 3.1, we need two more lemmas. 

Lemma 3.3. 

f a(x)SIv.SIvdx + f wrjdx = f aw(aiiVxi + CL22Vx2)dx 
2 Ot Jfi Jq Jq 

T / ((Jll^xi T 012^2)a;i (®12^a:i T 0,22Vx2^dx T / (tll2^a:i T 0-22W 2 )a:2 (®12'f^xi T (^22Xx2^dx. 
Jq J f2 

(3.21) 
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Proof. From (3.9) and (3.15) we derive 


—div{a{x)Vv)t = div{a{x){Vw + {xi — X 2 )e)) in'Dfil.). (3.22) 


Using (3.10)-(3.11), and taking into account (3.2), we obtain 

<—div{a{x){\/v)t,v>= / a{x)Vvt.Vvdx — / a{x)Vvt.iy vda{x) 

Jn Jon 

= j a{x)Vvt.Vvdx ~ y o,{x)Vv.Vvdx. (3.23) 

Similarly, we get by using (3.10), (3.17), and taking into account (3.2) 

< div{a{x)\/w),v > = — a{x)\/w.Vvdx + / a{x)Vw.u vda{x) 

.Jn Jdn 

= — a{x)\/w.Vvdx = — a{x)Vv.\/wdx 

J D. Jq, 

which can be written using (3.9) as 

< div{a{x)Vw),v > = — wrjdx. (3.24) 

Jn 

Similarly, using (3.10) and the fact that tti, U 2 satisfy (3.28)-(3.29), we get for fig = 
(0,L) X (eg, 77 - eg) 

< div{(ja - J^)a{x)e),v >= - iJa - X2)a{x)e.Vvdx 

J Qq 

= / + aw)a{x)e.Vvdx 

J Qn 


I no 

= / [div[a{x)Vv) + aw)a{x)e.Vvdx 

J Qq 

= / div{a{x)SIv))a{x)e.SIvdx + / awa{x)e.Vvdx 

J Qq '7 r^Q 

— I “1“ 0 ^ 12 '^X 2 )xii^l 2 '^xi “1“ (l 22 '^X 2 )^^ 

J Qq 

“t“ / {^^12'^xi H“ )ai2 ~\~ Cl22'^X2)dx 

J Qq 


+ 


/ aw{ai2Vx^ + a22Vx2)dx. 

J Qq 


Hence by combining (3.22)-(3.25), we get (3.21). 


(3.25) 

□ 


15 


Lemma 3.4. There exists a positive eonstant C sueh that 


— [ a{x)Vv.'\/vdx+ 2 f w{xi — X 2 )dx < C f a{x)Vv.Vvdx. 
Jq Jq Jq 


(3.26) 


The proof of Lemma 3.4 requires a lemma. 

Lemma 3.5. There exists a positive eonstant C sueh that 

rK-eo 


f °\Vv{0,X2)\‘^dx2 <C f \Vv\^dx 

t/ 6 q t/ 

f \Vv{L,X 2 )\‘^dx 2 < C [ \Vv\^dx 

J CQ J O 

f \Vv{xi,eQ)^dxi < C f \Vv\‘^dx 

Jo Jn 

f \Vv{xi, K — e())\‘^dxi < C ( \S/v'^dx. 

Jo Jn 


(3.27) 

(3.28) 

(3.29) 

(3.30) 


_ pK—eo 2 r® rJ^—eo 

Proof, i) Since V € C^(^o), we have / |VT(0, X 2 )Phx 2 = lim - / / iVvl'^dx 

Jen < 5^0 ^ Jo Jen 


5 pK-eo 


0 Jen 


It follows that we have for hi > 0 small enough 

pK—eo 2 /■'^l pK—^O 

/ |VT(0,X2)phx2 < ^ / / \S/v\^dx 

Jen Jo Jen 

which gives (3.27). 

In the same way we establish (3.28), (3.29), and (3.30). □ 

Proof of Lemma 3.4- We shall estimate the three integrals in the right hand side of 
(3.21). First, we obtain by applying Young’s ineqnality and using (1.2)-(1.3) 


/ aw{ai2Vxi + a22Vx2)dx < ^ / w‘^dx + ^ {ai2Vxi + a22Vx2y 

J Qq " J Qq J Qq 

[ w‘^dx + ^ I {aj^ + al^y^^'^dx 
Jn ^ Jn 


dx 


< 


< 


a 

2 


^ f w‘^dx + ^^^ [ a{x)Vv.Vvdx. (3.31) 
2 Jn 


m 
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Next, we have 


I ^12'^X2)xii^l2^xi ~\~ (l22'^X2^dx — / {Q'liyxi)xi{^12'^xi)dx 

J Qn J Qq 

/ {ailVxi)xAO‘22Vx2)dx + / {ai2Vx2)xA<^l2Vx^)dx + / {ai2Vx2)xA<^22Vx2)d- 

jQn JUn Jnn 


' f^o 
+ 


= Ii + I 2 + h + h- 


(3.32) 


Let us estimate the integrals /j. Expanding and integrating by parts, and using 
(1.2)-(1.3), (3.1), and (3.28), we obtain for a positive constant Ci 


h = 


I {an)x^au'^^dx + ^ [ auaui^^Jxidx 

J Hq ^ J Hq 

/ {aii)xiai2vl^dx-]- [ {aiiai2)xi^xidx 

J ^ d Hq 

2 j-K—eo 

{aiiai2vl^){L,X2)dx2 - - / (aiiai2U^J(0,X2)iix2 

^ Jen 


cK-eo 


< 


If ^ If ^ h? rK-eo 

- {aii)x,ai2vl^dx - - an{ai2)xivljx + — v^^{L,X2)dx2 

^ Jflo ^ Jflo ^ Jen 


^2 pK-eo 


^^{L,X 2 )dx 2 < Cl / a{x)Vv.Vvdx. (3.33) 

Jn 

Note that since u(0, X 2 ) = v{L, X 2 ) = 0 for 0 < X 2 < i^, we have 


< AN v^^dx + — 

JUn deg 


Vx 2 i^,X 2 ) = Vx 2 {L,X 2 ) = 0 for 0<X2<K. 


(3.34) 


Integrating by parts and using (3.34), we obtain 

I 2 = - [ aiiVxiia22Vx2)xidx = - I an{a22)xiVxiVx2dx - ]- [ aiia22{vljx2dx 

J r^o d s7q " J Hq 

= - [ aii{a22)xiVxiVx2dx+ 1- [ au{a22)x2'J^idx + ]- [ {aii)x2a22vl^dx 

JUn d Jq 2 Jq 


r*L 2 

{aiia22vi^){xi,K - eQ)dxi +- / 

/o d Jo 

Using Young’s inequality, (1.2)-(1.3), (3.1), and (3.29), we obtain from (3.35), for a 


1 1 _ 

-7; {ana 22 f^^){xi,K - eo)dxi + - {ana 22 vl^){xi,eo)dxi. (3.35) 

^ Jo d Jo 
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positive constant C 2 

h < AN f \vxj^\.\vx 2 \dx + AN f vl^dx + ^ [ vl^{xi,eo)dxi 
J £7n J £^0 " Jo 

A^ f 

\Nv\‘^dx-\ -/ ^{xi,eo)dxi < C 2 / a{x)Nv.Vvdx.{3.36) 

2 Jo Jn 


3AN 


Expanding and integrating by parts, and using Young’s inequality, (1.2)-(1.3), (3.1), 
and (3.30), we obtain for a positive constant C 3 


h = [ iai2)xiai2VxiVx2dx + I- [ o?2(^i)x2'^3: 

J S7q " J r2o 

= / {ai2)xiai2VxJ}x2dx [ {a\2)x2^xidx 

Jur, 2 


+ 0 / ia,i2Vxi)ixi,K-€o)dxi-- {a^2V ){xi,eo)dxi 


10 


< 


<Co 


- + vl^)dx - - ial 2 )x 2 vl^dx + — vl^{xi, K - eo)dxi 

J £^0 ^ J £^o ^ Jo 

/ a{x)VvN/vdx. (3.37) 

Jn 


Integrating by parts and using (3.34), and taking into account (1.2)-(1.3) and (3.1), 
we obtain for a positive constant C 4 


h = 


/ {ai2)xiO,22vi^dx + - / 012022(5^2)^1'^® 

J £lo J £^o 

= [ (oi2)xi0225^2*^® ~ o / (®12a22)a;i'yx2^® 

J £lo J £^o 

= 7 : (( 012 )xi 022 - 012 ( 022 )xi) 5 ^ 2 < Q / a{x)NvNJvdx. (3.38) 

J £7o J £^ 

To estimate the last integral in the right hand side of (3.21), we integrate by parts 
and use (3.29-(3.30) and (1.2). We obtain for a positive constant 

j {o-12'l^xi T O22'03;2)x2 (oi2'02;i T ^22'^X2)dx ~ X / ((oi2Xxi T <122^X2^ )x2dx 
J Qq ^ J Qq 

= -/ {ai 2 Vxi + a 22 Vx 2 f{xi,K - eo)dxi - - {auVxi + a 22 Vx 2 )‘^{xi,eo)dx 
2 Jo 2 Jo 
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< 


< 


1 ^ _ 

77/ {auvxi + a 22 Vx 2 )'^{xi,K - eo)dxi 
^ JO 

/ {{auVxif + ia22Vx2)‘^)(xi,K - eo)dxi 

Jo 

< / ((uxi)^ + {vx 2 )‘^)ixij ^ ~ ^o)dxi < Cs / a{x)Vv.'Vvdx. (3.39) 

Jo Jn 


Finally, combining (3.21), (3.31)-(3.33), and (3.35)-(3.39), we get for a positive con¬ 


□ 


stant C 

— [ a{x)Vv.Vvdx + 2 f w{xi — X 2 )dx<C f a(x)VF.Vi;(ix, which is (3.26). 
ot Jn Jn Jn 

Proof of Theorem 3.1. First, integrating (3.26) from 0 to f and letting h ^ 0, we get 

f a{x)\/v.\/vdx + 2 / f w{xi — X 2 )dxds < C f f a{x)Vv.Vvdxds. (3.40) 
Jn Jo Jn Jo Jn 

Next, we observe that since Ui G H{Xi) ™ Qy we have 

^'^(Xi - X 2 ) > 0 a.e. in Q. (3.41) 

Setting F{t) = ( f a{x)\/v.\/vdxds, we deduce from (3.40)-(3.41) that 
Jo Jn 

F'{t)<CF{t) VtG[0,T]. (3.42) 

Integrating (3.42), we get since F(0) = 0, 0 < F{t) < F{0)e^^ = 0 Vt G [0,r], or 

n a{x)Vv.Vvdxds = 0 Vt G [0,T]. 

! 

Using (1.2), we obtain Vu = 0 a.e. in Q. Taking into account that u = 0 on 
r 2 C dil and the connectedness of fl, we obtain u = 0 in Q. Going back to (3.3), we 
obtain rj = 0 a.e. in Q, which reads 


'»^ + Xi-X 2 = 0 in Q. (3.43) 

Multiplying (3.43) by w, we get aw'^ + w{xi ~ X 2 ) = 0 a.e. in Q. Taking into 
account (3.41), we obtain = 0 a.e. in Q, or ui = U 2 a.e. in Q. Finally, we obtain 
from (3.43) that xi = X 2 a.e. in Q. This achieves the proof. □ 


19 


Acknowledgments The second author is grateful to Prof. J. F. Rodrigues for kindly 
inviting him to the CMAF where he enjoyed excellent research conditions during the 
preparation of his Ph.D Thesis. 

References 

[1] J. Carrillo, On the uniqueness of the solution of the evolution dam problem, 
Nonlinear Analysis, Theory, Methods and Applications, 22(5), 573-607 (1994). 

[2] J. Carrillo and G. Gilardi : La vitesse de propagation dans le probleme de la 
digue. Ann. Fac. Sci. Toulouse Math. (5) 11, No. 3, 7-28 (1990). 

[3] J. Carrillo and A. Lyaghfouri : A hltration problem with nonlinear Darcy’s law 
and generalized boundary conditions. Advances in Differential Equations Vol. 
5, No. 4-6, 515-555 (2000). 

[4] Challal, S., Lyaghfouri, A. : A Filtration Problem through a Heterogeneous 
Porous Medium. Interfaces and Free Boundaries 6, 55-79 (2004). 

[5] S. Challal and A. Lyaghfouri : On a class of Free Boundary Problems of type 
div(a(X)Vu) = —div(H(X)x(u)). Differential and Integral Equations, Vol. 19, 
No. 5, 481-516 (2006). 

[6] S. Challal and A. Lyaghfouri : The Heterogeneous Dam problem with Leaky 
Boundary Condition. Communications in Pure and Applied Analysis. Vol. 10, 
No. 1, 93-125 (2011). 

[7] E. Dibenedetto and A. Friedam, Periodic behaviour for the evolutionary dam 
problem and related free boundary problems, Communs partial diff. Eqns 11, 
1297-1377 (1986). 


20 


[8] G. Gilardi : A new approach to evolution free boundary problems, Comm. part, 
diff. Eqns Vol. 4, 1099-1 123 (1979); Vol. 5, 983-984 (1980). 

[9] D. Gilbarg, N.S. Trudinger : Elliptic Partial Differential Equations of Second 
Order. Springer-Verlag 1983. 

[10] A. Lyaglifouri : The Inhomogeneous Dam Problem with Linear Darcy’s Law and 
Dirichlet Boundary Conditions. Mathematical Models and Methods in Applied 
Sciences 8(6), 1051-1077 (1996). 

[11] A. Lyaghfouri : The evolution dam problem with nonlinear Darcy’s law and 
Dirichlet boundary conditions. Portugaliae Mathematica 56 (1), 1-38 (1999). 

[12] A. Lyaghfouri : The evolution dam problem with nonlinear Darcy’s law and 
leaky boundary conditions. Ricerche di Matematica Vol. XLVII, Ease. 2, 297- 
357, (1998). 

[13] A. Lyaghfouri : A Regularity Result for a Heterogeneous Evolution Dam Prob¬ 
lem. Zeitschrift fur Analysis und ihre Anwendungen, Vol. 24, No. 1, 149-166 
(2005). 

[14] A. Lyaghfouri, E. Zaouche: IT-continuity of solutions to parabolic free boundary 
problems. Electronic Journal of Differential Equations, Vol. 2015, No. 184, pp. 
1-9 (2015). 

[15] A. Torelli, Existence and uniqueness of the solution of a non steady free bound¬ 
ary problem. Boll. U.M.L, 14-B(5), 423-466 (1977). 

[16] E. Zaouche: Existence of a Solution in a Class of Parabolic Free Boundary 
Problems. Submitted. 


21 


